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' Abstract 

Central production cross-sections of hard dijets with 2,1 or rapidity 
gap at Tevatron are analyzed in terms of diffractive ( "a la Good- Walker" ) 
'-^ ' and non diffractive fluctuations of the incident particles. The observed 

large factorization breaking and the unexpected high value of the 2 to 
1 gap cross-section ratio are explained in terms of scattering with and 
between the incident particles. 

1. Recent measurements of production of two large transverse momentum 
jets in the central rapidity region at Fermilab energies [1, 2] revealed a strong 
breaking of the Regge Factorization between the "no-gap" , "single-gap" and "two- 
gap" cross-sections. If the presence of the rapidity gap is interpreted as the colour 
singlet (c/. "Pomeron") exchange and its absence as the octet exchange, Regge 
factorization appears to be broken if 

= #2/i 7^ = #i/o > (1) 

where E , Ei and E 2 denote the cross-sections for "no-gap", "single-gap" and 
"two-gap" events, respectively, and Ri/j = Ej/Ej, by definition. 
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Experimentally [1], the ratio R2/1 (resp. R1/0) in (1) is estimated in practice 
by the ratio of dijet event rates per unit £ p (resp. £p), the fraction of momentum 
lost by the incident proton (resp. antiproton), measured as a function of XBj- 
Averaged over (small) XBj, the results are quoted to be R2/1 ~ 0.8 ± 0.26 for the 
L.H.S of (1) and R 1/0 ~ 0.15 ±0.02 for the R.H.S. [1]. The factorization breaking 
R1/0/R2/1 ~ 0.19 ± 0.07 [1] thus reaches a factor as small as 1 over 5. In itself, 
the high value of i? 2 /i = 0(1) is also a question to be understood. Our aim is to 
give a theoretical estimate of these ratios. 

In the present note we argue that the origin of these effects is basically the 
same as that observed in measurements of diffractive (virtual) photon-induded 
and hadron-induced processes [3, 2]. We show that it can be described in the 
same framework as done in Ref. [4] as due to the influence on dijet diffractive 
production of the soft scattering of the diffracted initial particles (protons and 
antiprotons at Tevatron). 

2. Let us recall briefly the argument given in [4] where the (single) diffractive 
production of two large transverse momentum jets was discussed. Following the 
idea of Good and Walker [5], one first expands the incident hadron state into a 
superposition of diffractive (eigen)states. 

\h >= ci|^i > +c 2 |^2 > +••• (2) 

A diffractive state is an eigenstate (by definition, in the subspace of states 
spanned by diffractive interactions) of absorption, i.e. we have 

T\$i >= > + Wilfa > (3) 

where T = 1 — S is the absorption operator and thus Aj is the absorption 
coefficient 1 of the state \ipi >. \<pi > describes the inelastic states which are 
at the origin of absorption. They are multiparticle states which do not show any 
rapidity gap 2 . Unitarity of the S-matrix implies 



|1 - Aj| 2 = 1 - \ui\ 2 = 1 - a nondiff , (4) 

where (Jnondiff is the non-diffractive cross-section. (4) is a special case of the 
generic relations for diffraction 

|1 — A| 2 = 1 — a ine l = 1 — CT nondi ff — o dd 
I — UJ 2 = 1 — (Jnondiff = 1 — Pinel + Odd (5) 

where Odd is the diffractive dissociation cross-section. Indeed, for a diffractive 
state, as seen from (3), a dd = 0. 

*At high energy Xi is dominantly real and positive. 

2 The first term in (3) is usually interpreted as "Pomeron exchange" where no colour is 
exchanged between \ipi > and the target. The second term would then describe all colour 
exchanges and thus \idi\ 2 represents the probability of this colour-exchange interaction. 
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Using (2) and (3), it is not difficult to express the diffractive transitions be- 
tween different hadronic states in terms of the absorption (A n ) and expansion 
(c n ) coefficients: 

<h'\t\h>=Y,^n{c' n Yc n . (6) 

n 

The physical content of the expansion (2) and of the formula (6) depends, of 
course, on the physical meaning one ascribes to the "diffractive" states >. 
Following [4] we could take them as states with a fixed number and transverse 
positions (in impact parameter space) of partons [6, 7]. In a modern language, 
we could equivalently consider QCD dipole states [8]. 

3. We are interested in three 3 processes: 

(i) p + p — > p + (C) + p 

(ii) p + p — ► p' + (C) + p 

(Hi) p + p -> p' + (C)+p' , (7) 

where (C) denotes a centrally produced system (at rapidity y) of two large trans- 
verse momentum jets and a "background" of the soft particles nearby, p', (p f ) 
represent the states which have no rapidity gap between the p (p) and (C) . Thus 
(i) corresponds to "two-gap" events, (ii) to "single-gap" events and (iii) to "no- 
gap" events, see Fig.l. 




Figure 1: Diffractive dijets with and without gaps. 



Now we have to write down the expansions of (7) into the diffractive states. 
Since we are interested in hard diffraction, i.e., in the process whose probability 
is small, we follow [4] and assume that the expansion is quasi-diagonal. Conse- 
quently we write 

3 It will also be useful to consider the process (ii)* symmetric to (ii) in the interchange of 
left and right moving projectiles, namely 

(m)* p + p ->■ p+ (C) +p' . 
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\p> = \g > + e\g + D > + e P \p + D > 
\p' + C> = - e*\g > + \g + D > + e'\p + D > 
\p + C> = - e* P \g>- e'*\g + D > + \p + D > (8) 

where e, ep and e' are small (and thus will be kept only up to first order) 4 . Here 
| g > denote a superposition of diffractive states representing a bunch of soft 
partons (close to their distribution in the left moving proton) and (D) is a su- 
perposition of partonic states consisting of a hard (i.e. small in transverse space) 
dipole at rapidity y and a number of soft partons with rapidities nearby (this sys- 
tem eventually decays predominantly into two large transverse momentum jets 
and the background of soft particles). 
Using (3) and (8) we thus obtain 



<p\t\p> = \p = \g 

<p' + C\t\ P > = e(\ { g + D) -\g) 

<p + C\t\p> = e P (X(p+D) ~ X g ), (9) 

where the transition matrix t is to be considered diagonal in both rapidity and 
impact parameter space. All quantities on the R.H.S. are to be understood as 
linear combinations of the corresponding quantities defined for truely diffractive 
parton states. 

Following [4] we now assume that \( 9 +d) an d \( p+ d) for diffractive states can 
be calculated assuming an independent scattering of their components [9], i.e. 

1 - \ 9+ d) = (1 - A 5 )(l - X D ); 1 - X {p+D) = (1 - Ap)(l - X D ). (10) 

When this is inserted into (9) we have 

< P ' + C\t\p>=€ X D (1-Xg) 

< p + C\t\p >= e P X D (1 - Xg) . (11) 

4. To obtain the cross-sections E lj2 for one-gap and two-gap events, see (1), we 
have to square the corresponding amplitudes, sum over the final hadronic states 
and integrate over the suitable rapidity and impact parameter intervals. In a 
first stage we shall discuss the differential cross sections a± t 2, before integration 
over rapidity and impact parameter variables. For given rapidity y and impact 
parameter position s, the summation extends for o~\ over all states (p 1 ) and (C) 
as the detailed final state is not measured. For 02 one sums only over all states 
of (C), because the final proton is identified. These differences imply different 

4 Thc relations between expansion coefficients follow from orthonormality of the states. 
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averaging procedures for the diffractive states. Taking this into account, we 
obtain 



<*i(s,y) = Hs,y)\ 2 \\ D (b-s,Y-y)\ 



<?2(s,y) 



\e P (s,y)\ 2 \\ D (b-s,Y-y)\'< 



\ g (b,Y) 



(12) 



—* 

where we have explicitely indicated the averaging procedures. Y, b are the rapidity 
and impact parameter distances between the incident particles. At this point one 
can observe that, as seen from (4), (5), 



1- 



u 2 (b,Y)} =l-uj 2 p (b,Y) = (l-\\ g (b,Y)} ) 2 +a dd . (13) 



It is important to remember that - since the averaging concerns only one vertex - 
add here corresponds to the single diffractive dissociation in one vertex. All this 
implies that 



a 1 (s,y)=[\e(s,y)\ 2 \\\ D (b-s,Y-y)\ 2 } [1 - co 2 p (b,Y)} 



(14) 



so that in the following one can omit the index av without running into confusion. 

One sees that the formula (12) for a 2 is not symmetric with respect to the 
interchange of the projectile and the target. To restore the symmetry we have to 
require that 



\\ D (b-s,Y-y)\ 2 = V P \ €p (b-s,Y-y)f 



(15) 



where, a priori, Vp is a vertex function of (s, y; b—s, Y^y) which must be symmetric 
with respect to exchange of the corresponding arguments. However, since the 
L.H.S. of (15) depends only on b— s and Y—y, Vp can depend neither on s nor on 
y. Symmetry implies that it must be a constant, depending only on the internal 
variables of the vertex. 

The formula (15) shows that the two-gap diffractive interaction can be equiv- 
alently understood as the elastic interaction of the projectile with a Good- Walker 
(Pomeron-like) fluctuation of the target. 

From (12) and (15) we deduce that 



<r 2 (s,y) = V P \ep(s,y)\ 2 \e P (b-s,Y-y)\ 2 \l-\ 9 (b,Y)\'< 
= Vp 1 \\ D (s,y)\ 2 \\ D (b-s,Y-y)\ 2 \l-\ g (b,Y) 



(16) 



5. We have to discuss now the no-gap events and thus to go beyong the Good- 
Walker argument, as the latter refers only to diffractive interactions. To this end 
we observe that the probability of a no-gap event to occur can be written as the 
product of (a) the probability of fluctuation of the projectile into \g' + D >, i.e. 
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|e| 2 , and (b) the probability of non-diffractive interaction of the central system 
(D) with the target, i.e. uj 2 d : 

a = \e(s,y)\ 2 cu 2 D (b-s,Y-y) (17) 

Since the same argument can also be used by exchanging the roles of the target 
and projectile we also have 

a = \e(b-s,Y-y)\ 2 u 2 D (s,y) (18) 

which implies that 

[u D (s,y)] 2 = VW,y)\ 2 (19) 

where V is another vertex function, symmetric with respect to exchange of the 
corresponding arguments. For the similar reasons as those following Eq.(15) we 
deduce that it is actually a constant. Thus we finally obtain 

a = V \ e (s,y)\ 2 \e(b-s,Y-y)\ 2 

= V- 1 u 2 D (s,y)u; 2 D (b-s,Y-y) . (20) 

A closer look shows that the vertex functions V and V P are identical. This 
can be seen by observing that the cross-section for one-gap events, <7i, given in 
(12), can be also calculated as a product of (a) the probability of fluctuation of 
the target into (p + D), i.e. |ep| 2 , (b) the probability of non-diffractive interaction 
of (.D) with the projectile, i.e. uj 2 d and (c) the probability that no non-diffractive 
interaction of the final antiproton with the projectile took place 5 , i.e. (1 — uOp) 
(c.f. (5)). Thus we can write (up = u p ) 

a 1 (s,y) = \e P (b-s,Y-y)\ 2 u 2 D (s, y) [l - u 2 p (Y)] . (21) 
Comparing (21) and (14) we have 

\e P (b-s,Y-y)\ 2 [cu D (s,y)} 2 = \e(s,y)\ 2 [X D (b-s,Y -y)] 2 . (22) 

This equation has a simple physical meaning. It says that the probability of 
a fluctuation of the projectile into the projectile+gap+D configuration is propor- 
tional to the elastic interaction in the projectile-D system, whereas the probability 
of fluctuation into no-gap configuration is proportional to the non-diffractive in- 
teraction in this system. It thus provides the interpretation, in the Good- Walker 
language, of the "Pomeron exchange" mechanism for gap creation. 

From (15), (19) and (22) one sees that indeed the vertex functions V and Vp 
are identical: 

V = V P . (23) 



5 Diffractive excitation of the projectile is allowed. 
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Using now (22) and (21) we obtain 6 

a^y) = V \e(s,y)\ 2 e P (b-s,Y-y)\ 2 [l-u; p (b,Y) 2 } 

= V- 1 \LU D (s,y)\ 2 \\ D (b-s,Y-y)\ 2 [l-u p (b,Y) 2 ] . 



(24) 



6. Putting together the formulae Eqs. (16), (20) and (24) for the differen- 
tial cross-sections, and even before evaluating directly the observable quantities 
mentionned in section 1, it is already possible to get important qualitative hints 
on the pattern of factorization breaking in diffractive dijet production at the 
Tevatron. 

Considering for instance, the differential 7 ratios r^/j = o-j/oj, we find 



r 2 /i 



n*/o 



\e P (s,y)\ 2 \l-X g (b,Y)\ 2 \X D (s,y)\ 2 \1 - X 9 (b,Y) 



Hs,y)\ 2 \ 

\tp(s,y)\ 2 
Hs>y)\ 2 



i-"|(y) 



I^d(s,2/)| 2 

\Ms,y)\ 2 



l-^JY) 



\uD{s,y)\- 



i-^(y) 



(25) 



where we make use of <r*, (cf. footnotes 3, 6) instead of <Ji in the second row for 
formal convenience, the qualitative conclusions being unchanged. 

Formulae (25) show explicitely that the factorization formula (1) is violated: 



vjt&y) _ 

<7o(s,y) <ri(s,y) 



<T 2 (S, y) [1 - CT n ondiff(b, Y)f 

X 



a inel (b,Y) 



(26) 



We have expressed the coefficients \ g , u p in terms of their physical interpreta- 
tion (4), (6). One sees that the factorization violating factor has the same origin 
as that which was shown in [4] to be reponsible for violation of factorization be- 
tween single diffractive processes at HERA and at FERMILAB. The numerical 
value may be, however, somewhat different, depending on the size of the cross- 
section for the soft diffractive dissociation in pp collisions (giving the difference 8 
between a ind and cr nondiff ). 

Another interesting remark individually concerns the ratios r 2 /i and T\ */ . 
Rewriting expression (25) as 



7-2/1 = 



\ e P\ 



1 Ginel 
&nondiff 



6 We can equivalently consider the symmetric a\ of o\ by interchange of left and right moving 
particles, namely 

01 (si 2/) = ai(b-s,Y-y) . 



7 Not to be confused with the integrated ones (1). 

8 In [4] this difference was neglected. If taken into account, the corresponding factor between 
the HERA and FERMILAB cross-sections should be 1 ^Vnondiff rather than 1 — tJi ne i- At small 
impact parameters, however, which are most important for numerical estimates, this difference 
is expected to be small. 
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|ep| 2 

ri */° = -J^T&y) i 1 - Vnondiff) ■ (27) 

In fact, we have to take into account that the specific factor |ep| 2 /|e| 2 is factor- 
izable since it reflects the probability ratio for a color singlet over an octet to be 
coupled to the projectile. Hence, we see that the essential factorization breaking 
factor is concentrated in the second ratio r±* /o and it has the same content as the 
corresponding factor between the HERA and FERMILAB cross-sections [4] . On 
contrary, the factorization breaking factor in r 2 /i is mild, since it depends only on 
the difference between 1 — Oi ne \ and 1 — <J non diff due to the diffraction dissociation 
contribution in the total cross-section. These conclusions will be confirmed by 
the following phenomenological application. 

7. To obtain numerical estimates of the cross-sections £j and of the measured 
ratios Ri/j, cf. section 1, one has to perform integration over the rapidity interval 

— * 

y and impact parameters b and s. 

For rapidity intervals, the integration is, in some sense, taken into account 
by the measurement itself, since event densities by unit of £ are given, hence, 
even if a full simulation would be welcome to match precisely the experimental 
conditions, the rapidity dependence is considered to be already integrated out. 

The impact parameter dependence has to be taken into account. This is 
difficult because several necessary elements are not known, so that -at best- only 
a rough estimate is possible. Only the elastic pp amplitude is known with some 
precision: 

\ p (b,Y)=a e- b2 ^; a = -pg- . (28) 

At v/i = 1800 GeV, a tot = 71.71 ± 2.02 mb and B el = 16.3 ± .5 GeV~ 2 [10], so 
that a 0.85. There is also the inequality (Pumplin bound [12]) 

\ 2 D < cu 2 D . (29) 

To simplify the discussion we assume tentatively that the shapes (in impact 
parameter space) of \ 2 D and uj 2 d are similar and take the Gaussian forms for easy 
integration: 

[X D (s,y)} 2 = r, u 2 D (s,y) = A exp[-s 2 /B D }. (30) 

where Bp is the slope in the elastic scattering of the central system (D). Guided 
by the results from hard diffraction of virtual photons [11], we estimate Bd to be 
in the region around 4 GeV -2 . 

Using (16), this allows to determine a 2 . The integration over d 2 s and d 2 b 
gives 
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where ( = Bp/ B e i. Using the obvious identity 

R 

R2/1 = R2/0R0/1 — -5 — (32) 
-rti/o 

one can thus calculate R2/1 from (31) and the measured ratio R\/q. 

The Fermilab measurements [1] give R1/0 ~ 0.15 ±0.02. Substituting all this 
into (32) and using the measured ratio R2/1 = 0.8 ± 0.26 we deduce that for 
( ~ 0.25, i] m 1 (with an error of about 30%). Although this value of rj is a 
perfectly acceptable one 9 , it should be noted that it saturates the unitarity limit. 
More precise data are needed, however, to consider this interesting property as 
really established. 

This result shows two points. First that our analysis can naturally explain a 
rather large value of the ratio R2/1 measured by the CDF coll. This seems to be 
a rather non-trivial result. As a sort of by-product we obtain the second point, 
namely the strong violation of the Regge factorization. 

To obtain the theoretical value of the ratio 

R = R2/1/ R1/0 — R2/0/ R-1/0, (33) 

one needs additional information about the soft single diffractive dissociation 
cross-section in one vertex (c.f. Eq. (4)). For a rough estimate we take 



< Aj(6) >= al e - b2/B ^s (34) 

which is a simplified version of the analysis by Miettinen and Pumplin [7] . From 
this and (28) we deduce the diffraction dissociation cross-section to be 

°dd = (« - l)<7 e j (35) 
where 10 k = B di ff/B el . This gives 

R ^ = " t 1 " 2 TTc + t4a) (36) 

9 The fact that we can take 77 ~ 0(1/2 — 1) is related to the fact that the central system 
D contains some soft partons apart from the hard dipole, e.g. inclusive diffractive production 
[13]. Therefore we would predict that in an experiment which will measure the "elastic" two 
jets e.g. exclusive diffractive production without accompanying soft hadronic radiation [14], 
one shall expect r)«l and falling with increasing transverse momentum of the jets. 

10 The ratio k > 1 does not imply that the slope in diffractive dissociation must be larger 
that in elastic scattering. This was explained already in Ref. [7]. The parameter Bdiff is 
responsible for the behaviour of the average of sum of the squares of the diffractive amplitudes 
(including elastic amplitude) in impact parameter space. From this information it is not possible 
to deduce directly the behaviour of diffractive dissociation channels in momentum transfer. As 
shown by Ref. [7], although the diffractive dissociation is more prcripheral in b space than 
elastic scattering, its slope in momentum transfer is smaller than that of elastic scattering. As 
explicitely derived and clarified in the parton picture of Ref. [7] , the slope in momentum transfer 
depends strongly on degrees of freedom which are other than the transverse ones and thus is 
not given by the fourier transform of (34) (cf. the comparison of Fig. 2 (for da/d 2 b) and Fig. 
5 (for da/dt) in Ref. [7]). 
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One sees that the ratio R is independent of 77 and thus completely determined by 
the value of k. To estimate k we observe that the total soft diffractive dissociation 
cross-section is approximately equal to the elastic one. Assuming the approximate 
factorization between the single-difffactive and double-diffractive cross-sections 
we can write for the sum of all contributions to the total soft diffractive cross- 
section 

o-difftot/o-ei = 2(k - 1) + (« - l) 2 = {k 2 - 1) w 1 

leading to k For ( = 0.25 this gives R xs 4, in agreement with the 

experimental value 5.3 ± 1.8. 

8. In conclusion, we have analyzed the factorization breaking observed in dijet 
diffractive production at Tevatron, using the same framework which allowed to 
explain the factorization breaking between HERA and Tevatron single diffractive 
processes. In both cases, the long distance diffractive interactions between the 
proton and the antiproton are the cause of the factorization breaking mechanisms. 

Interestingly enough, the pattern of factorization breaking in this framework is 
quite dependent of the number of rapidity gaps ( "no-gap" , "single-gap" and "two- 
gap"), in agreement with the experimental findings. The factorization breaking is 
strong for the comparison between "no-gap" and "single-gap" cross-sections, since 
it is related to the total pp absorption factor (4), while it is weak 11 for "single-gap" 
vs. "two-gap" ones, where it is related to the fraction of the inelastic diffraction 
to the total inelastic cross-section which is small. 

Thanks to the identification of a factorization breaking mechanism, an out- 
come of our approach is the numerical evaluation of the ratios Ri/j and in par- 
ticular R 2 /i which is found experimentally to be 0(1), that is surprisingly high. 
With quite reasonable values for the absorption parameters we find results in nice 
agreement with the data. Stimulating dijet production results, expected coming 
soon [16] from the Tevatron, Run II, will allow us to perform a more differential 
analysis of the mechanism we propose. 

Note added. While completing this study, the paper Ref. [17] has appeared, 
treating the same question within the rapidity gap formalism. We note an agree- 
ment between the two approaches. 
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